POINT MASS INSERTION ON THE REAL LINE AND 
NON-EXPONENTIAL PERTURBATION OF THE 
RECURSION COEFFICIENTS 



MANWAH LILIAN WONG 



Abstract. We present the construction of a probability measure 
dj with compact support on K such that adding a discrete pure 
point results in changes in the recursion coefficients without expo- 
nential decay. 



1. Introduction 



Suppose dfi is a probability measure on the unit circle <9B. We define 
an inner product and a norm on L 2 (<9B, d/i) respectively as follows: 



Using the inner product defined above, we can orthogonalize 1, z, z 2 , . . . 
to obtain the family of monic orthogonal polynomials associated with 
the measure d/i, namely, (& n (z, d/i)) n . We denote the normalized fam- 
ily as {ip n (z,dn)) n . 

Closely related to <& n (z) is the family of reversed polynomials, defined 
as ®n( z ) = zn ®n(l/z)- They obey the well-known Szego recursion 
relation 
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and a n is known as the n-th Verblunsky coefficient. The Szego re- 
cursion relations for the normalized family of orthogonal polynomials 
is 

p n +i(z) = (1 - \a n \ 2 )- 1/2 (zip n (z) - c^<p* n (z)) (1.4) 

These recursion relations will be useful later in this paper. 

Now we turn to a probability measure d'y on R. We can define an 
inner product and norm on L 2 (IR, d^y) as in (1.1) and (1.2), except that 
in this case it does not involve any conjugation. By the Gram-Schmidt 
process, we can orthogonalize 1, x, x 2 , . . . and form the family of monic 
orthogonal polynomials, (P n (x))^ =0 . Upon normalization, we obtain 
the family of orthonormal polynomials, (p n (x))™ =0 . These polynomials 
satisfy the following three-term recursion relation 

xp n (x) = a n+1 p n+1 (x) + b n+1 p n (x) + a n p n _i(x) (1.5) 

where a n and b n are real numbers with a n > 0. They are called the 
recursion coefficients of d*y. 

The main result of this paper is as follows: 

Theorem 1.1. There exists a purely absolutely continuous measure djo 
supported on [—2,2] with no eigenvalues outside of [—2,2], such that if 
we add a pure point x £ ^\[ — 2, 2] in the following manner 

dj{x) = (1 - P)d l0 {x) + /36 xo > (1.6) 

it will result in non- exponential perturbation of the recursion coeffi- 
cients a n (c?7o) andbni^djo). 

This example is of particular interest because of the following his- 
tory: back in 1946, Borg [TJ proved a well-known result concerning the 
Sturm-Liouville problem that in general, a single spectrum is insuffi- 
cient to determine the potential. Later, Gel'fand-Levitan [8] showed 
that in order to recover the potential one also needs the norming con- 
stants. 

Norming constants correspond to the weights of pure points and 
it is known that in the short range case (in orthogonal polynomials 
language, a n — l,b n — > fast), varying the norming constants will 
result in exponential change in the potential. 
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Moreover, when considering the effect of varying the weight of dis- 
crete point masses on orthogonal polynomials (both on R and 90), 
Simon proved that it will result in exponential perturbation of the re- 
cursion coefficients (see Corollary 24.4 and Corollary 24.3 of [15J). 

All the results mentioned above lend to a few the intuition that if the 
recursion coefficients a n — > 1 and b n — > fast, then adding a pure point 
will result in exponential change in the recursion coefficients. However, 
it turned out not to be the case! 

2. Tools Involved in the Proof 

2.1. The Szego Mapping. It turns out that one can relate measures 
supported on [—2, 2] with a certain class of measures on dD. 

Note that the map 9 \— > 2 cos# is a two-one map from <9B to [—2, 2]. 
Therefore, given a non-trivial probability measure d£ on <9B that is 
invariant under 9 — > —9, we can define a measure 

d 7 = Sz(d£) (2.1) 

using what is known as the Szego map, such that for g measurable on 
[-2,2], 

J g(2 cos 6)d£(6) = J g{x)d 1 (x) (2.2) 

Conversely, if we have a probability measure ft supported on [—2, 2], 
we can obtain a probability measure 

v = Sz" 1 ^) (2.3) 

on <90 by what is known as the Inverse Szego Mapping, such that for 
h(z) measurable on <9B, 

J h{9)dv{9) = j h (cos" 1 1) dry{x) (2.4) 

There are many interesting results about the Szego mapping (see 
Chapter 13 of [E]), but the only relevant one for this paper is the 
following by Geronimus jH] (see also Theorem 13.1.7 of [T4"]): 

Theorem 2.1 (Geronimus [10J). Let d^ be a probability measure on <9B 
which is invariant under 9 — >■ —9 and let dj = Sz(£). Let ot n = a n (d£), 
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a n = a n (dj) and b n = b n (dj). Then for n = 0, 1, 2, . . . , 

a 2 n+ i = (1 - a 2 „-i)(l - a 2 „) 2 (l + a 2n +i) (2.5) 

b n +l = (1 - «2n-l)a2n - (1 + a 2 n-l)«2n-2 (2.6) 

with the convention that a_i = —1. 

2.2. The Point Mass Formuula. We add a point mass ( = e lU) e <9D 
with weight < /3 < 1 to dfi in the following manner: 

dv=(l- P)dfi, + P8 u (2.7) 

Our goal is to investigate a n (di^). 

Point mass perturbation has a long history (see the Introduction of 
[17J). One of the classic results is the following theorem: 

Theorem 2.2 (Geronimus [TU1 HTj). Suppose the probability measure 
is defined as in (2. 7| ). Then 

^■^- (i-ffra,o (2 - 8) 

it 

ff„(*,c) = $3^<Pi(*) ( 2 - 9 ) 

3=0 

and a// objects without the label (du) are associated with the measure 
d/j,. 



Since $ n (0) = — a n _i, by putting z = into (2.8) one gets a formula 
relating the Verblunsky coefficients of d\x and dv. For more on the 
formula (2.8), the reader may refer to Nevai [T3l IT2] . and Cachafeiro- 
Marcellan [31 121 El H G3 . 

Using a totally different approach, Simon [T3] found the following 
formula for OPUC: 

a n {dv) =a n - g-^n+itC) a ^ ^$ + || ^ i(C)J ( 2 - 10 ) 

where q n = (1 - f3) + /3K n ((); a_i = -1. 
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Simon's result lays the foundation for the point mass formula. In 



[TBI HZ], Wong applied the Christoffel-Darboux formula to (2.10) and 



within a few steps from (2.10) proved the following formula for a n (dv): 

a n (du) = a n (dfx) + A n (C) (2.11) 

where 



An(C) 



:i-\an\ 2 ) i/2 <pn + i(cmo 



^n(C) = El^(0| 2 ( 2 - 12 ) 
3=0 



Formula (2.12) turns out to be very useful (see for example, [HI [T7J 

ESI). 



3. Outline of the Proof 

3.1. Case 1: xq > 2. We construct a measure d'jo with recursion 
coefficients (o n ) and (b n ) satisfying 



a n / 1 



b n = 
\ 2 = oo 



(3.1) 
(3.2) 



The measure d'jo is purely absolutely continuous and symmetrically 
supported on [—2,2], with no pure points outside [—2,2]. We scale it 
by a factor < y < 2 to form the measure dj y supported on [—y, y] C 
[—2, 2] (we will show the connection between y and Xq a bit later; see 



Then we use the Inverse Szego map on d'jy to obtain the measure 



dfiy. By looking at the Direct Geronimus Relations (2.5) and (2.6), we 



find necessary conditions for a n (djj, y ) so that both (3.1) and (3.2) hold. 

Since d^ y is supported on [— y, y] C [—2,2], we know that d\i y is 
supported on two identical bands. Besides, &\i y is symmetric along 
both the x-and y-axes because of the symmetry of d'jy and the Szego 
map. 

We add a pure point at z — 1 to d\i y to form the measure djx y and 
compute the perturbed Verblunsky coefficients a n {djx y ). 

Then we use the Szego map on djx y to obtain the probability measure 
dj y on K. Finally, we scale dj y to form the measure dj. 
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Note that if we have chosen y such that 

2 \x \ 



(3.3) 



then we have dq = (1 — f3)d r y Q + /3S Xo . 

As the final step, we show that for some constants C Xo ,D Xo (both 
dependent on x ) such that 

a n {dj) = a n (d l0 ) + ^ + o ( (3.4) 

bn(dl) = K(d l0 ) + ^ + o (J^j (3.5) 

3.2. Case 2: xq < —2. Everything in Case 1 will follow except that 
we add a point z = — 1 to dfi y instead. As we shall see later in the 
proof, dfiy is symmetric both along the x— and y— axes. Therefore, 
adding a pure point at z = —1 is the same as adding a pure point at 
z = 1 and then rotating the measure by an angle of 7r. 



For the convenience of the reader, here is a diagram showing all the 
measures involved. We will start from the measure d/u, y , and move 
along two directions: 

, scaling , Sz -1 , add z=l ,~ Sz" 1 ,~ scaling , 

"7o ^ — «7y ^ — "/-^ — "Ay — >* «7y — "7 (.<J-"J 
4. The Proof 

Let g?7o be a probability measure on M with recursion coefficients 



satisfying (3.1) and (3.2). 



This measure, supported on [—2,2], is purely absolutely continuous, 
and has no eigenvalues outside [—2, 2]. Moreover, if we write d^{x) = 
f(x)dx, f(x) is symmetric. 

Now we scale dj to form the measure d^ y defined by 

d ly {x) = d 1 {2xy- 1 ) 0<y<2 (4.1) 
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The measure dj y , supported on [—y, y] C [—2, 2], is purely absolutely 
continuous and the a.c. part of d^ y (x) is 

f y (x) = fC2.ru -) V y , y] (4.2) 

which is also symmetric. 

It is well known that scaling has the following effects on the recursion 
coefficients 

a n (dj y ) = (|j a n (d^ ) b n (dj y ) = (|j b n (d^) (4.3) 

Now we apply the inverse Szego map to d^y to form the probability 
measure dfi y on <9B, see figure below: 




Figure 1. Graph of supp (<i/i) 

The measure d[i y is supported on two arcs, [9 y , tt — 9 y ] and [tt + 
9 y , 2tt — 9 y ], with a.c. part 

w y {9) = w y {9) | [e y , v -0y] + u;„(0) | ^+e y ,2n-e y ] (4.4) 

where 

w„(0) = 2ir\sm(9)\f y (2 C os9)x [ e y , 7 r-e y] (9) (4.5) 

= cos- 1 (|) E (o, ^ (4.6) 

By Corollary 13.1.8 of [14], b n (d^ y ) = if and only if a2 n (dfi y ) = 0. 
Therefore, we can express the Verblunsky coefficients of dfi y as 

0,70,0,71,0,72,... (4.7) 



8 M.-W. L. WONG 

with Tj = tt2j+i- Moreover, by Theorem 13.1.7 of [Hj, we know that 

a2 n+\{dly) = (1 - «2n-l(^))(l - a 2 n(dfi y ) 2 )(l + a 2 n+l(d// y )) 
= (1-T n _i)(l+T n ) 

(4.8) 

Now we will choose a suitable family of r n G R such that the corre- 
sponding a n (d^ y ) satisfy both (3.1) and (4.8). 
Observe that by (4.8) above, 

a n+1 (d-f y ) 2 - a n (d^ y ) 2 = (1 - r„_i)(r n - r n _i) + (1 + r n _i)(r n _i - r n _ 2 ) 

(4.9) 

(4.10) 



Therefore, if we have an increasing family of r n < such that 
r n ^T 0O = -\ I - (-] <0 



then a n {d^ y ) /*■ y/2 and the corresponding measure &\i y . 
In particular, if we let 

then the goal is achieved. 

Next, we prove the following lemma: 



(4.11) 



Lemma 4.1. Let d\i y be the measure on dB) with Verblunsky coefficients 



as in (4.7) where for all large n, 

1 

T~n Too 7= 

m 



1 < r„o < 



(4.12) 



We add a pure point at z = 1 to d/j y as in (2.1) to form dfi y . Then for 
n = 2m or 2m + 1, A n (l) has the following expansion 



1 



A„(l) = - 
Therefore, 

a n (dfi y ) = 
where e m = o (m~ 3//2 ) . 







1 



2-Tno / m 3 / 2 



1 



m 



3/2 



(4.13) 






1 








f 2r 

^ 1 oo 




1 




m 3/2 



+ e m n = 2m 
+ e m n = 2m + 1 



(4.14) 



Proof. Since all the Verblunsky coefficients of d\i y are real, by induction 
on the recursion relation (1.4), 



n-l 



Vn(l) 




(4.15) 



3=0 



By (4.7), when n = 2m or 2m + 1, 

</?:(!) = ^(i) = 



m— 1 




(4.16) 



This formula will play a crucial role in the computation below. 

4.1. n is even. First, we compute A n (l) when n = 2m using the point 
mass formula (2.12). Let 

B n 

Then 



(1 - 7)7^ + «■„(!, 1) 



lim A 2m (l) = lim (1 - \a 2r . 



|2\l/2 



^2m 



,. ^2m 

lim — — 

m— >oo x5o 



(4.17) 
(4.18) 

(4.19) 



because «2 m = for all m. However, instead of computing this directly, 
we use the Stolz-Cesaro theorem (see [7]), which reads as follows 

Theorem 4.1 (Stolz-Cesaro [7j). Let (Tk)k, (®k)k be two sequences 
of numbers such that 6 n is positive, strictly increasing and tends to 
infinity. If the following limit exists, 



lim 

then it is equal to lim^oo Tk/Qk- 
First, note that — >• < 0. Thus, 

B m nK n {l,l) > |<^(l)| 2 ^oo 



(4.20) 



(4.21) 



by (4.16). Hence, it is legitimate for us to use Theorem 4.1 above. 

Let K n = K n {l,l) and y? n = ip n (l). Observe that (p 2m +i = f>2m- 
Therefore, by (4.16), 



B 



2(m+l) 



- B 2m - ¥2(m+l) + fir, 



2</>L 



(4.22) 
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A 2 (m+1) - A 2m - V J 2(m+l) ~ ¥ 2 m 



-2t„ 



l+T„ 



As a result, 



lim A 2m (l) = lim 



^2(m+l) — A 



2 m 



■m-^oo B 2 {m+1) — B 2m 

Next, we will prove that the rate of convergence is 



A 



2ml 



"Too + 



1 



+ 



m 



by computing the following limit 

lim m (A 2m (l) + T a 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



Recall the definition of A n (l) and the facts a 2m = and ^ 2m +i^2m = 
<p 2m - Thus, the left hand side of (4.26) can be expressed as X n /Y n 
where 



y 
Y 



+ T 00 K 2m ] 



K 2m -> oo 



We use the Stolz-Cesaro Theorem again. First, observe that 

1 Trn. . 2 



1 +T„ 



1 



l + r r 



(4.27) 
(4.28) 

(4.29) 



Then we compute 



X 2 (m+l) — X- 



2 m 



y 1 ^TTL 

vm + 1- 

1 + T m 



<^2m+ T oo 

s 



Vm + IK; 



2(m+l) 



2)» 



(1) 



(II) 



(4.30) 



Consider each term in (4.30) above. 

Vm + lj^ ~ V^+T(l - r m ) - v^(l + r„ 



(I) 



Y 2 (m+l) — ^2m 



Moreover, 



l+r„ 



(4.31) 



(II) = ^"1 + l(AT 2 (m+i) - K 2m ) + (v / m+ 1 - y/m)K 2m (4.32) 
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which implies that 
(II) 



Y 



m+l 



Y 

- 1 it. 



vm+l + (vm-j-l — \/m 



1 + r m K- 



m n 2m 



(4.33) 



Next, we show that lim^^oo K 2m /<^f m = — l/r^ by the Stolz-Cesaro 
Theorem. 

K2m 



lim 



lim 

\ m— >oo 


2 

^2(m+l) ~~ 


-^2(m+l) _ 


lim I 


1 — r 

1 ' m 


771— >OQ V 


1 + T m 






Too 





■^2m 
1 



(4.34) 



Combining (4.31), (4.33) and (4.34), we obtain 
lim m (A 2m (l) - (too)) = 1 

m— >oo 

Next, we are going to show that 



A 2m (l) = "Too + 



1 



+ o 



by computing the second-order term. We do so by proving that 

1 



lim m ( A 

m— >oo 



2m 







Let 
Then 

Pm+l ' 



m<p 2m + mTooK; 



p 



(m + 1) 



1 — 7V, 



m 



2m 



'film, 



mK- 



2m 



(4.35) 
(4.36) 

(4.37) 
(4.38) 



l 2m 



+ (m+l) ^2 (m+l) - K 2m ] + [(m + 1) - m] t od K 2 

- yjm + 1 [iC 2 ( m +i) - K 2m \ - (\/m + 1 - y/m)K 2m (4.39) 



Combining with previous results about Y m+ i — Y m and Ki m j^\ my we 
have 

Pm+l -Pn 



lim 







(4.40) 



which proves (4.36). 
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Next, we will obtain the third-order term by computing 

1 



L, = lim m 3/2 A 2 , 



^2m 



—7V, 



Let 



J m = m 3/ VL + rn 3/2 T O0 K 2m - mK- 



2m. 



By a similar argument as in (|4.34|) 

(m + 1) 



J'in+1 J in 



3/2 1 T m _ m 3/2 



1+7V, 



(4.41) 

(4.42) 



m 



'2m 



+ l) 3 / 2 ^ [K 2(m+1) - K 2m ] + [(m + 1) 3/2 - m 3 / 2 ] 

- (m + 1) [^ 2 (m+i) - K 2m ] - (m + 1 - m)Ar 2m (4.43) 
which implies that 

(4.44) 



L 3 = lim = 1 



Y 

- 1 rr. 



2r n 



4.2. when n is odd. We compute A n (l) when n = 2m + 1 using the 
point mass formula (2.12). Let A n and B n be defined as in (4.17) and 
d4~l8| . Then 

lim A 2m+1 (l) = (1 - M 2 ) 1 /* lim (4-45) 
We will use the Stolz-Cesaro Theorem again. Note that 



1 7Vn+l 1 I'm / 1 



and because v?2m+3 = <^2m+2, 

B 2 ( m +i)+i - B 2m+1 = 2ipl m+2 = 2 ( - - ipl m (4.47) 
Therefore, 

_ r (l _ \ T I2U/2 

lim A 2m+1 (l) = ^ ^ N = (4.48) 
"° V(l + r 00 )(l-r 00 ) 
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Next, we prove the rate of convergence by computing 

lim vM A 2m+i(l) + t-oo) = 1 (4.49) 

m— 5>oo 

Since a n G R, the recursion relation becomes 

(1 - \a n \ 2 ) 1/2 (p n+1 = <p n - a^<p* n = (1 - a n )(p n (4.50) 



Therefore, 



A m (1 ~ «2m+lV|m+l /, ^ x 
A 2m+ l(lJ = — = (1 - T m ) 



<P2m 



K 2m +l K 2m+ i 



(4.51) 



Let 

P m = Vm [(1 - T rn )(p\ m + r 00 K 2m+ i] (4.52) 

Q m = i^2m+l -> oo (4.53) 

Note that 

Qm+l — = -^2m+3 — ^2m+l = 2^2(m+l) (4.54) 



and 



Pm+i -P m = Vm + l(l - r m+1 )tpl {m+1) - y/m{\ - r m )ipl n 

s 



+ TooV^ + 1 [K 2m +3 ~ K 2m+1 ] + (Vm + 1 - ^/m)T 00 K 2m +i (4.55) 

V v ' 

(ii) 

Since (1 - T m )ip 2 2m = (1 + r ro )^ (ro+1) , 

(!) = Vm + 1(1 - r m+ i) - y/m(l + r m ) 

Qm+l 2 

Next, consider (II). We compute 

Um ^2(m+l) ~ ^2m = \l+r m 1 ) = - Too 



which implies 

n (II ^ n = -(1 + r 00 )(v / ^TT - v 7 ^) (4.58) 

Therefore, 

lim m (A 2m+1 (l) - (-Too)) = lim ^ = 1 (4.59) 
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Next, we will prove that 

A 1 ( X 

by showing 

Lf, = lim m ( A 2m+ i + 



As explained in ( |4.51[ ), it suffices to consider 

H m = m(l - + mr 00 K 2m+ i - s/mK 2m+1 



(4.60) 



(4.61) 



(4.62) 



H m+ x - H m = (m+ 1)(1 - r m+ i)^2( m+ i) - "7,(1 - T m )<y? 2 



(i) 



+ (m + l)r 00 ir 2m +3 - ^oo^2m+i -vWTTir 2m+3 + ^/mK 2m +i 



— \y — 

(II) 



(III) 



Since (1 - T m )(p 2 2m = (1 + TmV^m+i)' we have 

(I) _ (m + !)(!- r m+ i) - m(l + r m ) 

Qm+l Qm 2 



(4.63) 



(4.64) 



(II) 



roo(m + 1) + (m + 1 - m)r c 



-^2m+l 
' 2 ¥ , 2(m+l) 



(III) (-VmTT+ y/m)K 2m+3 



Qm+l Qr 



2 V 2 



2(m+l) 



This proves that L' = and thus (4.60) 
Next, we compute 



L' 3 = lim m 3/2 ( A 2m+ i 



(4.65) 



(4.66) 



(4.67) 



By similar arguments as in (4.64), (4.65) and (4.66), we conclude 
that 



L'=l + 



1 



2r n 



This concludes the proof of Lemma 4.1 



(4.68) 
□ 
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Finally, we apply the Szego map to this perturbed measure d\i y to 
form the perturbed measure dq y on [—2, 2], which is defined by 

%{x) = (1 - i)d ly (x) + 7^=2 (4.69) 

For the sake of convenience, we temporarily denote a n = a n {djx y ). 
Since b n (d^ y ) = 0, if suffices to consider 

It is more complicated with a n {d^ y ). Recall that 

a n +i{.d^ y ) 2 = (1 - a 2n -i)(l - a ln)( l + «2n+i) (4.71) 
and we know that 

a n +i(d ly ) 2 = (l-r n _i)(l + r n ) (4.72) 
Therefore, upon solving the algebra, we obtain 

a n+1 (rf%) 2 - a n+1 (rf 7 ,) 2 = 2(1 + ^ )m3/2 + ° (4-73) 
Upon scaling, we have 
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